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TECHNICAL MEMORANDUM X-64870 


THE DEFINITION AND SPECIFICATION OF THE 
NEAR-EARTH ENVIRONMENTAL CRITERIA 
FOR SPACECRAFT THERMAL DESIGN 

SUMMARY 


The variation of the earth’ s thermal and albedo radiation received by 
a near-earth orbiting space vehicle or space payload as a result of temporal 
variation of the earth- atmosphere is discussed. The time average effective 
value is introduced and defined as the average of the effective value of the 
albedo or thermal radiation affecting a thermal system in orbit. The time 
increment used determines the orbital segment over which the values are 
averaged. Time average effective values are the important values to the 
thermal designer. The values are independent of the spacecraft and depend- 
ent only on the character of the earth-atmosphere system. The time average 
effective values have a probability distribution which is dependent upon the 
time increment. A statistical study of current satellite data can give these 
probability distributions. With these distributions the thermal designer can 
define confidence levels on predicted temperature ranges which are com- 
patible with engineering models for use in design, failure probabilities, and 
spacecraft cost estimates. Use of the distributions in environmental criteria 
guidelines will improve current criteria in that the guidelines will be based 
on satellite data and will be compatible with engineering models. 


INTRODUCTION 


The problem treated herein centers on the values and definitions of 
the time-dependent thermal environmental design parameters for a space- 
craft in a near-earth orbit as employed by the thermal design engineer. 
Table 1 presents the current thermal environmental criteria parameters as 
given in Reference 1 for thermally emitted earth radiation, reflected albedo 
radiation, and the total direct solar radiation. The uncertainty involved in 
specifying the direct solar radiation impinging on a spacecraft is caused by 
experimental errors and not solar variation of energy output. However, the 



TABLE 1. PRESENT THERMAL ENVIRONMENTAL CRITERIA [ 1] a 


Variation of the Time Average Effective 

Values About Their Global Values 



Earth Emitted Thermal 



Radiation Variation 

Time Increment 

Albedo Variation 

(W/m 2 ) 

At < 0. 3 hr 

„ +0.30 

°' 30 -0.15 

237 +28 
237 -97 

0. 3 < At < 3 hr 

0.30 ± 0. 10 

+ 24 
237 -48 

At > 3 hr 

0. 30 ± 0. 05 

237 ± 21 


a. Solar constant = 1353 ± 13. 5 W/m 2 at 1 AU 


variation present in the radiation coming from the earth is a result of the numer- 
ous parameters which affect those values and which are time dependent. These 
time-dependent parameters, which encompass the geophysical and meteorolog- 
ical environments that lie below the orbit of the spacecraft, are the subject of 
this report. 

The present criteria have several defects which can, be eliminated 
through a judicial statistical study. The variation of the time average effec- 
tive values are now, at best, an educated guess. The lack of defined con- 
fidence levels for the values makes them incompatible with engineering 
models used in design, failure probabilities, and spacecraft cost estimates. 

The variations given are altitude independent and are not a smooth function 
of tlie time increment but are step functions. The albedo and thermal infrared 
radiation values are not interconnected. This is unrealistic and leads to ill- 
defined root-mean-square procedures for the thermal design engineer. As a 
starting point for the discussion of a procedure for improving the criteria, 
the fundamentals and nomenclature will be considered. 
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DEFINITIONS 


The concept for this report can be explained through the primary con- 
sideration of the variation of the time average effective value. The definition 
relating to this term will now be introduced f 2] . 


Time Average Effective Albedo 

If ds is an element of area of the earth- atmosphere 1 surface at P in 
Figure 1, the polar coordinates (0 O , 0 O ) are the zenith and azimuth angles of 
the incident radiation of wavelength \ upon the surface located at colatitude 
t] and longitude east C, , and (0, 0) are the zenith and azimuth angles of a 
reflected ray, then the spectral radiance 2 in the direction (0, <f>) is given at 
time t by [ 2] : 

N (£• rj> t, 0 O , 0 O , 9, 0) « j p (£, rf, t, 0 O , 0 O , 9, 0) cos 9 cos 9 0 H 

A, TT A A. 

(i) 

where H is the spectral irradiance 3 from the direction (0 O , 0 O ) in watt cm -2 

pm” 1 and p (£, rj, t, 0 O , 0 O » 9, <p) is the spectral bidirectional reflectance, 

i.e. , the fraction of the radiation incident from the direction (0 O , 0 O ) which is 
reflected in the direction (9, 0). 

The total reflected radiation per unit wavelength from the element of 
area ds is given by the integral over the hemisphere normal to the surface, 


1. A surface is assumed which has the properties of a combined earth surface 
and atmosphere. 

2. Spectral radiance is the radiant energy per unit of time per unit wavelength 
that leaves a surface per unit solid angle and per unit projected area of that 
surface [ 3] . 

3. Spectral irradiance is the radiant energy per unit time per unit wavelength 
that is incident upon a surface per unit area [ 3] . 
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Figure 1. Geometry of the earth-atmosphere system, the 
sun, and the satellite. 


4 


e. 


W^ds = 


2ir 7r / 2 

/ / N fe, r| , t, e 0 , 00, 0 , 0) Sin 0 d0 d0 

|_ 0=0 0=0 


ds . (2) 


The spectral directional reflectance r (0 O , 0 O , 4,17, t), defined as the frac- 
tion of radiation incident from direction (0 O , 0 O ) which is reflected, is 


<a , r - W* ds 

1 \( l 00 ’ ^ V ’ ^ " H COS0 O ds 


27T 7t/2 

= “ / / P, (£, Vt t, 0 O , 0 O , 0, 0) cos 0 sin 0 d0 d0 . 


7 r J \ 

0 = 0 0=0 


( 3 ) 


For a perfectly diffuse (Lambertian) reflectance, 


P x (t> V, t, 0 O , 0 O , 0, 0) = p (£, 77, t) , 


( 4 ) 


which is independent of the incident and reflected directions. From this the 
following relation is obtained 


r x (0 o , 0o> C» V , t) — p^(£>, 77, t) • 


( 5 ) 


Therefore, r is also independent of the direction of incidence and reflection. 

A 

The average or "total 1 1 directional reflectance over a range of wave- 
length and V 2 given by 
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^2 

J 0o, V> 0 

r (0 O> 0o, ?, U, t) = h 

X 2 

/ H x d x 

A, 

X 2 277 77/2 

Iff />\(t’V,t,6 o .<t>O’0'<t>)H x cos 0 sin 0 dd d</>dX 

- Ai 0-0 <M> (6) 

^2 

w / H x d x 

X, 

The value r is also referred to as the "local” albedo when the wavelength 
interval covers the spectrum and a Lambertian reflectance is assumed. 

The flux received by a spherical detector of radius a at an orbital 
height h above the earth of radius R is given by: 


<f ss , V t, h) = J ff N x [ e 0 (f, 17), <w, 17), <9cr, 1 ?), r, rj, t ] Afl R 2 sin f dX df dij ( 7 ) 

where R 2 sin f, df, d rj is the reflecting surface element, Afi is the solid 

angle subtended by the spherical detector at the element of area, and { £> , 

ss 

ru J is the subsatellite point Q, Explicitly, (0o> are time dependent since 
'S s 

they are measured from the subsolar point S which is defined by' the declination 
and Greenwich hour angle of the sun. The coordinates ( 8 , <p) and the solid 
angle A £2 = jra 2 ^a 2 can be found from the geocentric coordinates of the satel- 
lite and of the surface (£,, 17 ) . The distance yo. is from ds to the spherical 
detector and is given by 

/x - \/(R + h ) 2 + R z - 2 R(R + h) cos O' 1 

where 0* is defined in Figure 1. 

Following the assumption that the earth has a homogeneous earth- 
atmosphere surface of constant Lambertian reflectance r , the flux received 
would be 


$ = f J f A £2 cos 9 cos 0 O H^R 2 sin t, df, dr] d\ 


( 8 ) 
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Since 


does not depend on 4 or r ) , a new set of geocentric coordi- 


nates (0’ , 0 f ) can be used to specify the surface element depending only on the 
relative position of P and Q . Therefore, 


<T‘ = 


= P 


“ II 

7T J J 


cos 0 cos 0 n R 2 sin 0 r tra 2 d0’d0 t 
(R + h) 2 + R 2 - 2R(R + h) cos 9 ' 


(9) 


The total reflectivity has been used and is defined by 


K\ ax 

p = 

/«* dX 


The solar constant S has also been introduced. 


00 

S = f H dX . 
0 A 


It is seen from Figure 1 that 


cos 9 a = cos 9 ' cos 6 + sin 9 * sin 9 cos 0 ’ 

s s 


and 


cos 9 


(R + h) cos 9 ' - R 


Hence, 


9 ' 


max 77 


A . r * (cos O’ cos 6 + sin d* sin 6 cos 0 ) 

V* = 2a 2 SpR 2 j j * s - 3 ^ f(R + h) cos 0 - R] sin 6' d$' d0 f , 

n J n [R 2 + (R + h ) 2 - 2R(R + h) cos 6'] 


(IQ) 
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assuming the area of the earth below the detector is completely sunlit without 
the terminator visible where 

0' S 1 - 9 
max 2 s 


6 f is defined by the 0 T which specifies the surface elements on the local 
max 

horizon for the detector cos 0 ' = R/(R + h) . The integral in equation ( 10) 

max 7 

is the albedo form factor for a sunlit earth. 4 

The diffuse Lambertian reflectance for the earth is the albedo used in 
thermal design. Hence, for a specific time t , the effective albedo for an 
orbiting detector would be 


PJt) 


J J f P’ '■ ®o. 0) cos 0 cos 00 sin f df dTJ^- dX 

O' 

JTiaX TT f rr 

r r (cos# cos# + sin# sin# cos$)[(R + h)cos# - R] sin# d# d 0 

2S / / , , r* T: : 

Q 0 [R 2 + (R + h) 2 - 2R(R + h) cos # j 


(n) 


This defines the albedo as used for the models at time t . But as the detector 
moves in its orbit, the effective albedo changes because of changes in (l) the 
surface reflectance, (2) the area seen by the detector, and (3) the relative 
position of the sun, surface, elements, and the detector. Hence, the term 
'’time average effective albedo” (TAEA) is introduced and defined by 


t+At 


P (f* At) — J p ( t’) dt 1 , 

8 t 8 


4. The following approximation can be used for the albedo form factor [4l 


max 7T 


* s i 

0 0 


(cos # cos # + sin # f sin # cos (j)) 


[R 2 + (R + h) 2 + R(R + h)cos #'] 


[(R + h) cos B* - R] R 2 sin#' d#' d <t> 


= 2 i 


-J± !(tff 


1 + il 

R 


1 ^0.84 + 0.16 exp ^-60 — cos0 s ^ 
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It is the effective albedo seen by a detector taken over a time increment At . 
The statistical variation of this quantity can be used as the working criteria to 
determine the variation of the radiation due to albedo. 


Time Average Effective Thermal Radiation 

If ds is an element of area of the earth- atmosphere surface at point P 
and (8 , <p) are the zenith and azimuthal angles of an emitted ray of wavelength 
X , then the spectral radiance of the thermal emitted radiation at colatitude r\ 
and a longitude t, is given at time t by [ 5] 

rn BB 

(tj, £, t, 8, 0) = e x (9, 0, r ] , £,, t) N a [T(tj, l, t)] watts/m 2 sr 

( 12 ) 

where e(8, 0, tj, t) is the directional emissivity of the surface element ds . 
The blackbody spectral radiance is given by 

\ BB - 2c2frx " s [ exp (ot) - x ] ' (13) 

where c is the speed of light, T is the temperature of the body* k is 
Boltzmann’s constant, and h is Planck's constant. 

The energy radiated in all directions per unit area is given by the 
radiant emittance 


2ir ir/2 

W m (r?, I, t) = f J N “(r), t,, t, 8, <p) cos 9 sin 9 68 6<p 


(14) 


0 0 


XR 

If the surface is Lambertian, then N is independent of 9 and 0 ; hence, 

A 


t) = *N A m (r], i, t) . 


(15) 


Assuming a blackbody earth-atmosphere system in the infrared is equivalent 

to assuming that the system is Lambertian in that spectral region and € = 1 . 

A 
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The total blackbody energy radiated per unit area over the entire spectral 
range is given by 


W(n, t) = / W dX = / 


27TC 2 ff dX 


0 xS [ ex p(is) - 1 


= qT 4 


(16) 


where a = 27r 6 /l5'K 3 k 4 c 2 = 5. 729 x io _8 W/m 2 °K 4 . Hence, at time t and loca- 
tion (tj, £,) , a temperature T defines the energy radiated. The energy radiated 
in the direction 0 per unit solid angle per unit area is given by 


jUt = qT 4 
7T 


(17) 


The thermal energy received by a spherical detector of radius a at an orbital 
height h above the earth of radius R is given by 

,i , t) = ff cos 9 A ft R 2 sin £ dt, dr? ( 1 8) 

SS SS 7T 

where the integral is over the area of the earth-atmosphere surface seen by the 
detector and Aft (= na?/ja 2 ) is the solid angle subtended by the spherical de- 
tector from the element of area ds. The subsatellite earth coordinates are again 

(tj , £ ) . Using the geometry discussion from the albedo section, 

s s s s 


*IR (t? t £ f t ) = 2a 2 R 2 f f 2 aT 4 (e', 0!)[(h + R) cos £ - R] sin ffdO’d*’ 
ss ’ ss ’ 0 0 [R 2 + (R + h) 2 - 2R(R + h) cos 0 T ] 1 

(19) 

If T is independent of (0', </>’), then 


$ IR (t) = 2ia 2 trT < 


1 


n/ 2Rh + h 2 

(R + h) 


( 20 ) 


10 



For a satellite, the observed or effective temperature T is determined at a 
specific time t ; hence, 


* IR (h . t, h) = * 
ss ss 


IR 


(t) 


( 21 ) 


The effective blackbody radiation per unit area is given by 


IR(t) = <r T 4 


27ra 2 


^ IR (t) 

2Rh + h 2 
" (R + h) 


( 22 ) 


IR(t) defines the effective thermal emitted radiation for a model of the earth- 
atmosphere system at time t . As the detector moves in its orbit, the effec- 
tive thermal emitted radiation changes as a result of the surface temperature 
and its characteristics, i.e. , changes of atmospheric conditions, topography, 
season, and time of day. Hence, the term "time average effective thermal 
emitted radiation" is introduced and defined by 


t+ At 

IR = / IR(t') dt T . (23) 

t 


It is the effective radiation_seen by a detector over a time increment At . The 
probability distribution of IR as a function of At will be used to define the 
criteria in order to estimate limits on the temporal variations of IR . 


PARAMETERS AFFECTING THE TIME AVERAGE EFFECTIVE VALUES 

Geophysical Parameters 

The various parameters that influence the time average effective values 
fall into three categories: geophysical, meteorological, and geometrical. The 
geophysical parameters are those that are independent of time and refer to the 
characteristics which arise from the physical scattering and reflection prop- 
erties and the emissivities of the earth’s surface and atmosphere. A detailed 
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discussion of these quantities is beyond the range of this report; however, a 
general summary for albedo is given in Table 2 [6] and Figure 2. The albedo 
of clouds covers a range from 0.10 to 0.80, land features cover a range from 
0. 05 to 0.45 except for snow- and ice-covered areas, and water covers a range 
of 0. 02 to 0.20. The bidirectional reflectance from the earth surfaces and 
clouds could be obtained. However, the scattering of a Rayleigh atmosphere 
with aerosols presents a different picture. Snoddy concluded, "Due to the 
many variables and unknown factors, it is doubtful that a complete theoretical 
model for the reflected and scattered radiation from a planet and its atmosphere 
will be formulated in the near future" [7] . The effect of a Rayleigh atmosphere 
without aerosols can be inferred from the fact that above the subsolar point a 
homogeneous earth with p = 0 will have an effective albedo due to a Rayleigh 
atmosphere of about 0.08 [7] . Meteorological effects would drastically change 
this simple geophysical model. 

The geophysical factors that influence the emitted thermal radiation 
vary, depending on surface conditions. The emission properties are a function 
of the surface temperature, e.g. , desert, ocean, and polar region. At the 
atmospheric transmission region of 11 pm the following brightness temperatures 
were obtained from Nimbus 4 [8] : 


Sahara Desert 

321°K (102°F) 

Mid-Latitude Ocean 

280°K (44°F) 

Antarctic 

198°K ( -103°F) 


This illustrates the temperature variation over the earth. 

The water content in the air changes the thermal emission of the 
atmosphere due to the absorption bands (e.g. , H 2 0) . The various absorption 
and emission bands in the atmosphere complicate modeling by a simple geo- 
physical atmosphere in the thermal region. As an example of the brightness 
temperature of the atmosphere, the CO 2 band at 15 pm emits at 215° K (-76°F) . 
At such a temperature, the atmosphere is warmer than the antarctic surface 
temperature given previously. The total effective temperature is dependent 
on both the atmosphere and surface. 


Meteorological Parameters 

The meteorological parameters affecting the time average effective 
albedo and thermal emitted radiation are a function of time and location on the 
earth. Because the nature of the atmosphere is continuously changing, weather 
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TABLE 2. SUMMARY OF REFLECTANCE DATA FOR EARTH SURFACE 

FEATURES AND CLOUDS [6] 


Reflecting 

Surface 

Magnitude and Other Spectral 
Characteristics 

Angular Distribution of 
Reflectance 

Total Reflectance 

Soils and Rocks 

Increases to 1 pm 
Decreases above 2 pm 

Backseattering and forward 
scattering 

Sand has large forward scattering 
Loam has small forward 
scattering 

5 to 45 percent 
Moisture decreases 
reflectance by 5 to 
20 percent 

Smooth surfaces have 
higher reflectance 
Diurnal variation 
Maximum reflectance 
for small sun angles 

V egetation 

Small below 0.5 pm 
A small maximum bump at 0. 5 to 
0. 55 pm 

Chlorophyll absorption at 0.68 pm 
Sharp increase at 0.7 pm 
Decrease above 2 pm 
Depends on growing season 

Backseattering 

Small forward scattering 

5 to 25 percent 
Diurnal effects 
Maximum reflectance 
for small angles 
Marked annual 
variation 

Water Basins 

Maximum at 0. 5 to 0.7 jum 
Depends on turbidity and waves 

Large back and forward 
scattering 

5 to 20 percent 
Diurnal variation 
Maximum for small 
sun angles depends 
on turbidity and 
waves 















TABLE 2. (Concluded) 


Reflecting 

Surface 

Magnitude and Other Spectral 
Characteristics 

Angular Distribution of 
Reflectance 

Total Reflectance 

Snow and Ice 

Decreases slightly with increasing 
wavelength 

Large variability depends on 
purity, wetness, and physical 
condition 

Diffuse component plus mirror 
component ' 

Mirror component increases with 
increasing angle of incidence 

Variable 25 to 80 per- 
cent 

84 percent in Antarctic 
74 percent in Ross Sea 
ice 

30 to 40 percent in 
White Sea ice 

Clouds 

Constant from 0. 2 nm to about 
0. 8 pm 

Decreases with wavelength above 
0.8 ftm, showing water vapor 
absorption bands 

Pronounced forward scattering 
with small backscattering 
Minimum for scattering angles 
of 80 to 120 deg 
Fogbow for scattering angle of 
143 deg 

10 to 80 percent 
Varies with cloud 
type, cloud thickness, 
and type of under- 
lying surface 
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Figure 2. Local albedo variation [9, 10 ] . 


is highly variable and depends upon the temporal and spatial position over the 
earth. At any particular time the average cloudiness for the whole earth is 
between 50 and 60 percent. 

The local albedo should experience a diurnal variation as a result of the 
diurnal variation of temperature and sun angle for a particular location. The 
local albedo increases as the solar zenith angle decreases. The diurnal varia- 
tion of the albedo observed with Tiros 4 low resolution sensors is shown in 
Figure 3 t S ] . 

From climatological consideration, seasonal variations of albedo are 
to be expected. Figure 4 presents the seasonal variations of albedo as deter- 
mined from the 0.55 to 0.75 nm data from Tiros 7 [8] . Both sets of data were 
adjusted to give a mean quasi-global albedo of 32. 2 The seasonal variation is 
±0.05 from the global average. 

The latitudinal variation of the albedo is shown in Figure 5 [ 111 . These 
data were obtained from OSO-II by treating thermal control coatings as radiom- 
eters. The values were obtained during the period March 4-8, 1964. All 643 
data points lie between 33°N and 33°S latitude and are evenly distributed over 
all longitudes with the exception of the region between 40° W and 80° W. The 
OSO-I data and London' s (1957) latitudinal averages are also shown. Though 
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ALBEDO 










the data are limited, they do indicate the global latitudinal variation resulting 
in part from the decreasing solar elevation angle, snow and ice cover near the 
poles, and increased cloud cover associated with large-scale weather activity. 

A probability distribution function Z of the albedo for these data is also plotted 
in Figure 5. The distribution covers a selected temporal and spatial region for 
the earth. A smooth curve through the histogram indicates the probability, 
but, because of the difficulty of reading the data from the latitudinal figure, 
this curve only represents the trend. This is exemplified by the fact that the 
authors obtain an average albedo of 0.26, while the investigators report an 
average of 0.28. The albedo measurements range from 0. 10 to 0,62. The 
global average was 0. 29 when corrected. 

A representative variation of the albedo as seen by a detector is shown 
in Figure 6. These data of April 11, 196 5, were obtained by Linton [12] using 
thermal data from Pegasus I. The first half of the orbit shown has no cloud 
formation associated with it, and the albedo recorded is that of the surface 
features and the effects of a Rayleigh atmosphere. The central peak at longi- 
tude 80 deg corresponds to a cloud mass. 

The data from April 11, 1965; May 8, 1965; and March 20, 1965, were 
combined to form the probability distribution, also shown in Figure 6. This 
random sampling does agree in the gross feature with Figure 5 however bad 
the statistical sampling. 

The spatial variation over the globe is shown in Figure 7 for a 15-day 
period and was obtained by averaging the Nimbus II data over a grid field of 
5 by 5 deg of longitude and latitude at the equator and 5 by 2. 5 deg at 60 deg 
latitude [13]. The albedo of the earth system between 40°S and 40°N during 
the period May 16-31, 1966, is shown for a partial swath around the earth. 

An albedo distribution for the swath around the earth was generated as shown. 
Raschke [ 13] presents the reflected solar radiation and the emitted infrared 
radiation over the entire globe from the Nimbus II meteorological satellite 
during the period May 16 to July 28, 1966. In general, the albedo of the polar 
cap is in the range from 0.60 to 0. 90. Toward the equator the albedo decreases 
rapidly with decreasing cloudiness, showing minima of less than 0.20. The 
Sahara and Arabian deserts have albedo of more than 0.30 to 0.35. The global 
albedo was from 0.291 to 0.306. 

An Apollo 8 photograph of the Western Hemisphere shows the striking 
variation of albedo over geographic regions (Fig. 8) . Salient global cloud 
features can be obtained from Reference 14, which is a photographic summary 
of the earth T s cloud cover for the year 1967. 
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Figure 7. Frequency distribution of the albedo from 
Nimbus II between +40° N and -40° S during 
the period May 16-31, 1966 f 13] . 



Figure 8. Apollo 8 photograph of the Western Hemisphere. 


The diurnal variation of the thermal emission tends to peak in the early 
morning and decrease through the remainder of the day [ 6] . The diurnal varia- 
tion on a global average is ±28 W/m 2 . For a temperature change from 273° K 
( 32° F) to 288® K (59®F), the emission increases 76 W/m 2 from 318 W/m 2 . 
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The mean annual and seasonal variation of the emitted infrared radiation 
of the earth- atmosphere system observed from the first generation meteoro- 
logical satellites is given by Vonder Haar [6] : 

Thermal Emitted Radiation (W/m 2 ) 

Northern Southern 
Global Hemisphere Hemisphere 


Dec. -Feb. 

223 

223 

230 

Mar. -May 

230 

230 

223 

June-Aug. 

230 

237 

223 

Sept. -Nov. 

237 

237 

237 

Annual 

230 

230 

230 


The variation with latitude of the thermal radiation is shown in Figure 9 
for a time span from June through August. These data were obtained from 
Tiros VII in 1963. The thermal emitted radiation is lower at the poles. The 
variation of the absolute magnitude at most latitudes is very small. The largest 
seasonal differences occur in the polar regions and are about ±70 W/m 2 f 6} . 

To illustrate the variation of the thermal radiation seen by a detector, 
OSO-I and OSO-II data were used to derive a representative distribution (Fig. 

IC ) . The peak frequency is 330 W/m 2 with sideband widths of +70 and -90. 

Since the thermal radiation is correlated to the albedo through the heat 
balance equation, a general correlation is expected between the effective values. 
- r0ni Tiros VII data a correlation between the 8- to 12-^tm channel and the 0. 55- 
to 0. 75-fxm channel has been obtained: 


A albedo channel 
A IR channel 


0.0039 , 


where A indicates an incremental change in value. The data shown in Figure 
11 were obtained by using the longitudinal maps of Bandeen [ 2] . The data given 
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Figure 9. Tiros VII thermal and albedo radiation latitude 
data for June- August 1963 [ 2] . 


in Figure 7 were transformed into a distribution of thermal emitted radiation 
using the albedo and IR relationship and the following empirical formula: 

Thermal Emitted Radiation = 237 - 0.0039 (albedo -0.308) 


The results of this transformation are shown in Figure 12. The distribution 
about the mean, 237 W/ m 2 , is shown. 

The spectral distribution for the reflected albedo and thermal emitted 
radiation is shown in Figure 13. The meteorological effects change the spectral 
characteristics; however, the bounds indicated in the figure are representative 
of the fluctuation. 


Geometrical Parameters 

The geometrical parameters are derived from the relative locations of 
the sun, earth, and detector. These give the height of the detector, the angles 
for the bidirectional reflectance, and the various solid angles. 
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TCHmeHATURE, 








Figure 11. Tiros VII albedo-thermal radiation correlation. 


The area of the earth seen from a height h is given by 
j 4 ttR 2 = 2tt( 1 - cos 0 T ) R 2 ~ 2 tt ^ 2 R 2 = 7r(0'R) 2 (24) 

where Afi is the solid angle of the area as viewed from the center of the 

earth. If 0 is the angle from the normal to the local horizon where 9 ' = 
max 

0* (Fig. 14), then 
max 



0 


max 


2L 

2 


0 ’ 

max 


— 

2 


cos 


R + h 


(25) 
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EMPIRICAL EQUATION 


IR = 237-0.0039 (A-0.308) 



Figure 12. Empirical transformation of Nimbus II albedo data to 
obtain a thermal distribution function. 


For a 43 5 -km orbit, 9 =19.1 deg and /3 =70.9 deg. The effect of 

max max 

varying /3 keeping h constant is examined in Table 3 and Figure 15. For a 
Lambertian earth, halving the solid angle of the field of view should half the 
radiation received by a detector. If the solid angle is halved and the above 
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Figure 13. Spectral radiance functions for the albedo and thermal 
radiation emerging from the earth-atmosphere surface. 
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Figure 15. Effects of 0’ (/3, h) on surface area. 



approximation is used, then it is seen that half of the radiation (/3 A/2) is 

max 

coming from an area with a radius of 565 km about the subsatellite point. How- 
ever, an area of radius 2445 km influences the total effective albedo. 

To approximate the maximum rate of change the effective albedo can 
obtain through geometry, one assumes that the albedo p changes from 0 to 1 

G 

as the detector at height h moves 29 ' (Fig. 16). 

max 

Since the period of the orbit is 


P 


2tt 


(R + h) 3 


GM 


(26) 
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where cos = R/R + h. For h= 435 km (235 n.mi.), P - 93.6 and 

dp/dt ^ 1.7 x 10~ 3 /sec. 

Increasing the height will decrease the rate by increasing the time for 

the detector to move 20 1 — ► tt „ The area below the detector also 

max h— 

becomes larger and the average effective albedo is influenced by a larger area. 
But, as h increases, the influence of the albedo compared to the direct solar 
thermal input decreases. Figure 17 presents the view factor, defined by 

* received^ pS7ra2 for a s P here of projected area 7ra 2 and a Lambertian earth 
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Figure 17. The geometric view factor ($Ara 2 Sp) for reflected 
Lambertian radiation incident on a sphere versus altitude 
for a solar zenith angle of zero [15] . 


as a function of height [see equation (10) ] . Similarly, the view factor for the 

IR 

thermal radiation is defined as $ . ./lira 2 where I is the energy 

received t t 

radiated per unit area at the earth-atmosphere surface. 
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At 10 4 km the albedo input even for p = 1 is less than 0. 10. The 
value 10 4 km for h will be used to define the limit height for near-earth 
orbits. 


DERIVATIONS OF THE TIME AVERAGE EFFECTIVE ALBEDO 


The time average effective albedo depends on the geophysical, meteoro- 
logical, and geometrical inputs for the temporal and spatial regions under con- 
sideration. Therefore, the TAEA, 

At 

p g (t, At) = / P e (t’) dt' , (28) 


fluctuates between certain limits depending on At , the orbit, and the earth- 
atmosphere system in the swath of the detector. The distribution of the time 
average effective albedo for given parameters (such as the orbital height) 
versus the time increment is needed in thermal design. Several procedures 
for determining the TAEA will be discussed. The mean, lcr , 2ff , and 3a 
values (a is the standard deviation) are the quantities needed by the thermal 
design engineer because these quantities characterize the probability distribu- 
tion. The la , 2a , and 3a symbols indicate that the proportion of the values 
found within the interval of the mean plus and minus na (n= 1, 2, 3) are 
68. 3, 95.4 and 99. 7 percent, r espectively. This does not mean that the prob- 
ability distribution of TAEA Z [p (t, At)] is a standard normal distribution 
function. e 


Theoretical Distribution 

The following is a theoretical derivation of the probability function for 

the time average effective albedo. Physically, the TAEA cannot have values 

greater than 1. 0 or less than 0. The physical distribution is required to have 

a 0 probability at 0 and 1. The restricted range requirement does not allow a 

normal distribution function to approximate Z(p ) directly. Hence, a trans- 

e 

formation is introduced. The transformation of p^(At) to x by 
P (At) - 0.30 

x - -■ 6 ' (29) 

rp [i - p ( m) i 

v e 
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gives values at p (At) = 0.30 of x = 0 . 0 , at p (At) 1 an x value of +®, 
e 0 

and at p (At) ^ 0 an x value of . Hence, one can reconcile the normal 

distribution approximation using the parameter x . For the normal distribution, 
the probability of finding an albedo such that a ^ x b is given by 

P(a < x < b) = f e“ (x-p) 2 /2a 2 dx ( 30 ) 

a \ ! 2 ir o 

where a 2 is the variance (and a is the standard deviation), 


+ CO 


<x ! - / (x - „)■ 


e - (*-(*) Vac 2 a* , 


2 It (7 


ix is the mean, 


+ oo 


1 (x-u) 2 / 2 a 2 

-<*> n / 2 ir or 


p - / X — e dx , 


and 


P ( — co < X — co ) = 1 


The integrand of P(a^x^b) is the probability function Z(x) . Figure 18 
shows d ie normal distribution function Z and the corresponding values of 
P (At) under the transformation. In the above case, p = 0 and cr = 1 . Using 

the inverse transformation 


P e (A0 


x — 1 + \/(x — l) 2 + 1, 2x 
2x 


(31) 


the confiden ce value s in Table 4 ar e obtain ed. The probability distribution 

function Z [p (At)] in terms of p (At) is shown in Figure 19 along with the 
e e 

position of the la , 2 a , and 3 a values. 
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Figure 18. Normal (gaussian) probability function showing the 
corresponding p" locations under the transformation. 


TABLE 4. CONFIDENCE VALUES IN TERMS OF x AND o (At) 

e 


X 

(Minimum) 

— 
p (Maximum) 

G 

Probability (%) 

1 

0. 163 

0. 548 

68. 3 

2 

0.108 

0.711 

95.4 

3 

0.080 

0.793 

99. 8 


Because of the similarity of the OSO-II and Pegasus I data in Figures 5 
and 6 to the Z(p ) just derived, this function of Z[p (At)] shall be assumed 

G 0 

for the moment to describe the distribution of p^(At) at At= 0.- Further- 
more, it is assumed that at a time increment At= 2 hr , the function Z (p^) 

is normal with 3cr = 0.10, and as At-*-® it becomes a normal function with 
3cr = 0.05 . A 3cr of 0. 05 is similar to actual seasonal variation. The points 
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Figure 19. Probability distribution function for the effective 
albedo at an altitude h 0 and a time increment At. 


are used to curve fit Z [p (At)] to the function g 0 ± g t exp f ~k(At) ] for 

l(j, 2(j, 3cr values as a function of At. The results, with the present criteria 
values for comparison, are shown in Figure 20. The approach is not neces- 
sarily optimum, but it is a possible method of derivation of the appropriate 
criteria values, and the curves do indicate the physical behavior expected with 
a minimum effort. 


Satellite Derivation 

In an attempt to define the albedo probability distribution from satellite 
data, a simple procedure was developed to analyze Applications Technology 
Satellite (ATS) III photographs. An outline of the procedure is as follows: 

1. The geographic surface features were assigned albedo values and a 
2. 0 by 2. 0 deg grid was created. 

2. Visual inspection of an ATS photograph led to assignment of cloud 
albedo in grid form. 
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Figure 20. The 1, 2, and 3 cr confidence levels for maximum and 
minimum time average effective albedo. 



3. The two grids were merged. 

4. A Lambertian earth is assumed and, hence, the sun angle was 
neglected. 

5. Assuming a height and inclination, the position of the orbiting detec- 
tor was calculated as a function of time. 

6. The time average effective albedos were calculated. 

7. Finally, the lcr , 2cr , 3cr values were determined. 

A similar procedure for steps 1 through 4 was discussed by Bartman 
[ 2] ; however, neither the TAEA nor the usefulness of the results to the thermal 
design engineer was discussed. 

Figure 21 presents a portion of the cloud matrix between ±60 deg latitude 
and 0 to 70 deg longitude. The values of the albedo of the clouds were assumed 
according to the type and visual indications. Figure 22 presents the surface 
features where water was assumed to have an albedo of 0.10 and land one of 
0. 30. Several orbit paths are shown. The assumptions for the albedo of water 
and for underestimating the cloud cover lead to an accumulation of 0. 10 values 
in the albedo distribution for the earth-atmosphere system. This distribution 
is shown in Figure 23. 

For eight orbits and 400 points on the earth-atmosphere surface. 

Figure 24 shows the relative frequency of the calculated TAEA for At= 0. 5, 

5. 3, 16. 0, and 42. 5. The orbit was assumed to have a 50-deg inclination and 
a 556-km altitude. The resulting lcr and 3cr limits are given in Figure 25. 
The roughness of the calculation did not allow the actual 2cr and 3cr values 
to be determined. The 3cr given is the 100 percent probability limit. Figure 
26 shows the same results except that the inclination is 18 deg. Because of the 
lower albedo values near the equator, the mean global average is lower. The 
general trend is, again, the same as from the transform results. The error 
bar on Figures 25 and 26 indicates the probable error of the data points given 
the assumed albedo matrices for the surface and cloud features. This is an 
example of the type of data that could be extracted from satellite data. 


Satellite Data Study 

The general procedure to determine the thermal design parameters from 
satellite data is as follows: Convert existing satellite data from selected orbits 
for which both instantaneous values of the earth thermal and albedo radiation 
can be obtained into data which can be easily handled. Use these data to 
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Figure 22. Surface grid for the earth. 
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Figure 23. Albedo distribution for the earth-atmosphere surface. 


calculate the instantaneous values which are summed over various time incre- 
ments of consecutive parts of an orbit to obtain the time average effective value 
(TAEV). Determine the associated confidence levels by a statistical study of 
the different TAEV T s. Analyze the accumulated data for the effects of orbit 
inclination, angle between orbital plane and earth-sun line, time of year, and 
orbital position. Evaluate the results and present them in a form that can be 
used for design purposes. 

The second generation meteorological spacecrafts of the Nimbus series, 
specifically n and HI, seem to provide sufficient data to deduce information on 
the time average effective values. The Nimbus spacecrafts carried radiometers 
which were always pointed directly toward the earth. The orbits were nearly 
circular at about 1100 km and were almost polar. The angle of inclination of 
approximately 80 deg provided an orbit which remained oriented along the 
earth-sun line ( i. e. , sun-synchronous) . 
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Figure 24. TAEA distribution for At= 42. 5, 16. 0, 5. 3, and 0. 5 minutes. 
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Figure 25. Time average effective albedo distribution for 
50 deg inclination (8 orbit sum). 


The radiometers and spectrometers onboard provided data on the 
thermally emitted earth radiation and the reflected solar radiation emerging 
from the top of the atmosphere. The experiments on board Nimbus II and III 
which are of interest are listed in Table 5. The two channels of particular 
interest on Nimbus II are the medium resolution infrared radiometer (MRIR) 

5 to 30 jum and 0. 2 to 4 /^m channels. The 5 to 30 jum channel provides the 
radiance of 95 percent of the earth thermal radiation, while the 0. 2 to 4 jiim 
channel provides the radiance values of 99 percent of the solar reflected radia- 
tion. Because the data of the MRIR on Nimbus II did not provide complete full- 
earth coverage and the instrument operated only 3 months, the desined thermal 
design parameters cannot be completely recovered from the data. The data of 
Nimbus III should be employed for better statistical results. 
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Figure 26. Time average effective albedo distribution for 18 deg inclination (12 orbit sum). 



TABLE 5. NIMBUS SPACECRAFT SYSTEM 



Nimbus II 

Nimbus HI 

Activities 

May 15, 1966, launch 

May 22, 1966, operational 

July 26, 1966, 15 percent reduction of 

HRIR a data . 

July 29, 1966, MRIR data lost 

Nov. 15, 1966, HRIR data lost 

Apr. 14, 1969, launch 
Apr. 22, 1969, operational 

July 22, 1969, IRIS C data lost 

Aug. 12, 1969, noisy HRIR data 

Aug. 28, 1969, MRIR° data reduced 
Oct. 1-15, 1969, full coverage of earth 
by MRIR 

Jan. 31, 1970, MRIR pictures ter- 
minated 

Sept. 25, 1970, altitude control failure 

Orbital Data 

Height — 1179 by 1095 km orbit 
Inclination — 79. 689 deg retrograde 
orbit 

Period — 108 min 
Sun- synchronous 

Height — 1101 km 

Inclination — 80. 1 deg retrograde 

Period — 107 min 

Sun-synchronous 

Experiments of Interest 

1. MRIR Spectral Bands (pm) (50 km 
resolution) 

6. 4-6. 9 atmosphere H 2 0 and 
cirrus cloud mapping 
10-11 surface and cloud top 
temperatures 

1. MRIR Spectral Bands (pm) (50 km 
resolution) 

6. 5-7.0 atmospheric H 2 0 and 
cirrus cloud mapping 
20-23 atmospheric H£> and 
cirrus cloud mapping 



MU- 


TABLE 5. (Concluded) 



Nimbus II 

Nimbus III 

Experiments of Interest 

14-16 stratosphere temperatures 

10-11 surface and cloud top 

(Concluded) 

5-30 total thermal radiation and 

temperatures 


heat budget 

14.5-15. 5 stratospheric 


0. 2-4 daytime cloud mapping, 

temperatures 


albedw and heat budget 
2. HRIR Spectral Bands (gm) ( 10 km 

0. 2-4.0 daytime cloud mapping, 
albedo, and heat budget 


resolution) 

2. HRIR Spectral Band (g m) ( 10 km 


3. 5-4. 1 nighttime cloud coverage 

resolution) 

3.4-4. 2 nighttime surface and cloud 
cloud top temperatures and 
cloud mapping 

0. 7-1. 3 daytime cloud mapping 

3. IRIS Spectral Band 5-20 ^tm ( 150 
km resolution) 

Atmospheric temperature profile, 
O a , water vapor, and surface 
temperature 


a. HEIR — High resolution infrared radiometer 

b. MRIR — Medium resolution infrared radiometer 

c. IRIS — Infrared Interferometer Spectrometer 



Nimbus III did not have the 5 to 30 fim channel which provides direct 
information on the total earth-emitted thermal radiation. However, one can 
correlate the HRIR data of Nimbus II and III and the IRIS data of Nimbus III 
with the data from the MRIR' s of both spacecraft to provide data for a procedure 
to obtain satisfactory thermal data from Nimbus III. The Nimbus El data pro- 
vide a longer time of observation. 

A program to determine thermal environment parameters from satellite 
data is under study. The program wEl employ the Nimbus II and Nimbus III 
spacecraft experiment data obtained by the MRIR, the HRIR, and the IRIS 
experiments. 

Computer programs using the Nimbus data in digital format are to be 
developed to allow a statistical study of the data to determine the following 
near-earth thermal design parameter objectives for the earth’ s thermal and 
reflected solar radiation: 

1. Objective 1. The probability distribution function of the time average 
effective values and the 1, 2, and 3cr confidence limits (i. e. , the 68. 3, 95.4, 
and 99. 7 percent confidence level limits) for various orbital time integration 
intervals At’s. 

2. Objective 2. Global average value for the data maps processed under 
objective 1. The global average value for each map generated under objective 1. 

3. Objective 3. The effect of the variation of the orbit height and the 
orbital time integration interval and their interrelationship on the TAEV sigma 
values. The study will be performed by calculating the TAEV sigma value 
limits for various orbital heights . 

4. Objective 4. The relation of the earth-emitted thermal radiation 
and the albedo values and their sigma level confidence limits for simultaneous 
data will be determined by interrelating the 5 to 30 and 0. 2 to 4 pm channels of 
the MRIR on Nimbus II to allow the determination of the procedures to calculate 
the total thermal inputs into a spacecraft of both albedo and earth thermal radia- 
tion. The interdependence will be presented to allow the design engineer to 
calculate the root- me an- square averages of the heat input and specify these 
procedures. 

The interdependence will be in a form such that given a time average 
effective albedo for At , both the high and low values of the expected earth 
thermal radiation can be obtained for 1, 2, and 3<r confidence levels r i. e. , for 
a time average effective albedo ( At) , the earth thermal radiation has a 

ncr confidence of being below b^+ b^ exp (-b^At)] . 
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5. Objective 5. The distribution of the rate of change of the TAEV f s 
will be determined for the various orbital time integration limits. 

6. Objective 6. The relation of the sun-earth-spacecraft position 
(phase relation) on the time average effective albedo values will be determined 
for phase angles from 0 to 90 deg. This will give the dependence of the direc- 
tional reflectance of the albedo on the solar zenith angle. 

7. Objective 7. The effect of inclination on the sigma value confidence 
limits will be determined. 

The Nimbus II MRIR and HRIR data will be analyzed for various numbers 
of orbits, and the Nimbus TII MRIR, HRIR, and IRIS data will be analyzed to 
provide the correlation between the four thermal radiation channels of 6. 5 to 
7.0, 20 to 23, 10 to 11, and 14. 5 to 15, 5 of Nimbus III MRIR and HRIR channels. 
The correlation will provide a procedure to determine the total earth thermal 
radiation from Nimbus III data. 

Using the above correlation, the Nimbus n and III MRIR data will be 
transformed by integrating the data over areas of 160 km by 160 km into grid 
maps for the earth-emitted thermal radiation and albedo radiation levels. A 
grid map is a matrix of elements for which the location in the matrix is given 
by earth latitude and longitude and the value of the matrix element is the average 
value of the Nimbus data over an earth -atmosphere area of approximately 160 
km by 160 km. Each grid map is representative of the values obtained during 
a 24-hour period over the entire earth. Conversion for geometry and spectral 
band approximations are to be given. The earth is assumed to be a Lambert 
reflector for the solar radiation. Each grid map will be representative of the 
Nimbus data obtained for a 21 -hour period covering the entire earth. The days 
considered will be separated in time as equally as possible over a typical year 
or portions allowed by the lifetime of the experiments. For portions of the 
earth not covered by the Nimbus data during a specific day, appropriate inter- 
polation values will be estimated using data obtained on consecutive days or 
adjacent areas. These grid maps will model the earth albedo and thermal 
radiation for the various days. 

The constructed grid maps (i. e. , a matrix of albedo and thermal values 
of the earth- atmosphere system) will be employed to calculate the time average 
effective values of both the albedo and earth thermal radiation. The orbital 
node angle will be randomly distributed to average out any effects which it might 
introduce. The TAEV T s will be calculated using a set of orbits covering the 
entire modeled earth per grid map. The second half of the consecutive orbits 
should be displaced by one-half of the longitude displacement for one orbit. 
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For At = 0 the TAEV will be calculated at 30-sec intervals for determin- 
ing the other TAEV values and the effects of height, inclination, and orbital motion. 

The sigma level confidence limits will be calculated from the distribu- 
tion of the TAEV. The upper and lower confidence level limits of 68. 3, 95.4, 
and 99.7 percent for each orbital time integration interval given in objective 1 
will be determined. 

The results of the confidence level limits study will be presented in a 
tabular form giving the earth thei*mal radiation and albedo values at which the 
three confidence level upper and lower limits lie for each of the 100 days and 
the values for the average of the results for 3-month periods, A parameterized 

function of exp[ -c/ U> ^ At] , where C lt C 2 , and C 3 are param- 

eters and (u, l) specify the upper or lower limit parameter, will be fitted to 
each confidence level for the values of At. The display of the result will be 
in a form similar to Figure 20. 

The results of the global average per grid map and the total global 
average will be tabulated for the earth T s thermal and albedo values obtained 
from the study. 

The effect of the orbital height on the effective time integration limit 
will be given in a parameterized form so that, given the sigma level, confidence 
limits for any height can be determined from the results of a specific height. 

The parameters can be given in the form: 


- A 

= TTt 

a ‘ " TTt { Zlp e (dt)l 

These parameters will allow an 
increase in the time increment. 


At 


At 



increased height 


in orbit to be related to an 
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Since the bidirectional effects have been averaged out and the earth is 
assumed to be a Lambertian reflector, the effect of solar position is to be 
included. The limb darkening effects for thermal radiation will also be included. 


The rate at which p^ (At) and IR(At) can change with respect to time 

is important in creating various theoretical heat load inputs which are allowed 
under the given parameters. These distributions will be presented in the 
results. 

The relation between the thermal and albedo values will be included. 

This is a result of the fact that cloud cover increases the albedo but decreases 
i.he JR received because the radiating layer seen is higher in the atmosphere. 
This will affect die root-mean-square procedures employed by the designer. 


The specular ocean reflections will be included as singularities at the 
appropriate conditions. 


APPLICATION PROCEDURES 


In the past, the thermal design engineer has dealt mainly with slow 
response systems and systems which could be overdesigned, but now there are 
many thermal design applications in which the tolerances cause the uncertainty 
in the thermal evaluation to approach significant magnitude. For large-scale 
environmental control systems in spacecrafts, a 5 to 8° C (10 to 15° F) uncer- 
tainty in the albedo and IR effects has significant impact. 


Figure 27 shows the temperature of a shell versus orbital position with 
the surface density as a parameter. The thickness Ax necessary for a sur- 
face density p is given by: 

S 


p s = 


p V 

V 


~ 7rr J 

O 


- -ir(r - Ax) 3 
4 Trr 2 


(32) 


Table 6 gives the thickness for a 1-meter radius shell for aluminum and gold. 

It is seen that a small surface density corresponds to thin elements such as 
detectors and solar cell elements. As shown in Figure 27, as the mass of the 
surface increases, the temperature excursion in orbit decreases f 15] . Also, 
as the surface density increases, the calculated mean temperature increases to 
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10.0 LB/FT 2 
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40 60 120 180 200 240 280 320 360 

ORBITAL POSITION (DEGREES) 

Figure 27. Temperature history of an orbiting shell [15]. 


TABLE 6. SHELL THICKNESS AND DENSITIES 


p s 

(lb/ ft. 2 ) 

p s 

( kg/ m 2 ) 

10.0 

237,0 

3.0 

79.0 

1.5 

39.5 

0.75 

19.8 

0.25 

6.58 


Ax - A1 
= 2700 kg 


36. 7 cm 
12.2 
6.1 
3.0 


Ax — Ag 
= 19 300 k 


5. 1 cm 
.7 
.85 
.43 
14 












the limiting temperature which is given by 


T 

mean limit 


Table 7 gives the mean, maximum, and minimum temperatures shown in 
Figure 27 [15]. 


4 /rp 4 rp 4 

/ max min 

J 2 


( 33 ) 


TABLE 7. MINIMUM AND MAXIMUM TEMPERATURES FOR SHELL 


p s 

1 

Mean 

Max 

Min 

t/(T<. + 
' min 

T 4 )/ 2 

max 

(kg/m 2 ) 

[°C (° F) 1 

r°c ( °f) ] 

r°c (° f)j 

r°c c 

’F)] 

6. 58 

6.7 (44) 

67. 8 (154) 

-71. 1 (-96) 

22. 0 

(72) 

19. 8 

16.7 (62) 

59.4 (139) 

-35.0 (-31) 

23.3 

(74) 

39. 5 

19.4 (67) 

47. 8 (118) 

-10. 0 (14) 

23. 1 

(74) 

79. 0 

21.7 (71) 

37.8 (100) 

6. 7 (44) 

23.5 

(74) 

237. 0 

22. 2 (72) 

28. 9 (84) 

19.4 (67) 

24.3 

(76) 


The reaction of a system to the environment depends on the environment 
and on the design of the system. The albedo and IR criteria should be given 
so that they are independent of the system but are capable of being combined 
with the parameters of the system to give its thermal history. The significance 
of the time average effective albedo can be seen in this context by introducing 
a function called the "thermal time constant of the system." The temperature 
is assumed to be of the form 


T = T 0 e" tA (34) 

where r is the thermal time constant [ 16] . If it is near equilibrium, the 
thermal time constant can be defined as 
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T 


(35) 


me 

E 

16A(rc m T 0 * 


then 


r = p 400 sec m 2 /kg . (36) 

s 

Hence, for the shells considered earlier, 


r = 2fr, 3 hr 

p - 237 kg/ 
s 

= 8. 78 

- 79 

= 4. 39 

= 39.5 

= 2.11 

= 19.2 

= 0. 7311 

= 6.58 


For the lowest surface densities, the values of r are approximately equal to 
the time necessary for a half revolution. Hence, the time average effective 
albedo for At =* 30 min would be of particular importance. For a shell with the 
following parameters: 

e = 0. 0 
i = 33 deg 

h = 250 n. mi. (463 km) 

p = 1, 5 lb/ft 2 (39. 5 kg/m 2 ) 
s 

S = 1396 W/m 2 (443 BTU/hr/ft 2 ) 
a = 1.0 
e = 1. 0 
a /e = 1. 0 , 


the minimum and maximum values of the temperature for various albedos are 
given in Table 8. 
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TABLE 8. MINIMUM AND MAXIMUM TEMPERATURES FOR 

VARIOUS ALBEDOS 


Albedo 

Maximum Temperature 

r°c ( ° f ) j 

Minimum Temperature 
f°C (°F)| 

0. 20 

29.4 (85) 

-35.6 (-32) 

0.35 

33.9 (93) 

-38.9 (-38) 

0. 50 

37. 8 (100) 

-42. 2 (-44) 

0.80 

47. 8 (118) 

-49.4 (-57) 


The probability function for the time average effective albedo would indi- 
cate the chance of seeing albedos such as those in Table 8, or, equivalently, 
the chance of having those maximum and minimum temperatures. For a partic- 
ular daytime pass, the instantaneous average albedo can be any value between 
the limits given by the probability function of TAEA. 


A question that arises is, how long can a specific value p (t, At= 0) * 

6 

be seen assuming a certain Kcr value (K = 1, 2, 3)? The period of time is 
just the value of At which gives 


p (t, At) = p (t, At = 0)* . (37) 

G G 


For example, if K = 3 and p^ (t, At= 90) = 0.35, then one could have 

P g ( t, At= 45) = 0.474 for 45 min and p^ (t, At= 45) = 0. 226 for 45 min as 

long as these values are within the 3cr limit s for the distri bution. Further- 
more, the value during 15-min intervals for p (t, At= 15) could flip-flop 

G 

between 0 . 508 and 0. 192. This example represents a hot case for 3 <j since 
P (t, °° ) = 0. 35. Figure 28 represents various ways that such a hot case for 

a pass could be introduced into a thermal design program for thermal analysis. 
Each segment is such that if extended it would give an 0. 35 average over a 
period of time. 
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It, At-OI 



The instantaneous change is not necessary. A finite slope up to a maxi- 
mum of 1.6 x 10~ 3 sec could have been introduced for a 463-km orbit. 

The general limitations in parameterized form for the time history of 
the albedo are summarized in Figure 29. In this figure an example of a time 
history is given which considers those limitations. In generating fcriteria, 
various ’’hot" and ’’cold'’ examples of time histories can be generated as guide- 
lines for the thermal designer. The albedo, though, must be correlated with the 
thermal radiation. The parameter of the time average albedo could be displayed 
as in Table 9. Alternately, since for a specified cr limit 


p ( t, At) 


Z=K<r 


P e (t, At = «) + 


, „ -b Kl) At -c. (u ' 0 h/R 

a. u ’* e 1 e 1 

1 


(38) 


where (u, 1) are the upper and lower limits and a and a 11 have opposite 
signs, the values for a. , b , c. could be tabulated for season, inclination, and 

sigma limits. Similar results would be obtained for the thermal value. Then 
the thermal designer, knowing the system and expected orbit, could find the 
albedo and IR values and range that would be encountered by the system. 


CONCLUSION 


In summary, the variation of the earth thermal and albedo radiation 
received by a near-earth orbiting space vehicle or space payload as a result 
of temporal variations of earth-atmosphere has been discussed. The term 
"time average effective value" was introduced and defined as the average of the 
effective value of the albedo or thermal radiation affecting a thermal system 
in orbit over a given orbital segment defined by an increment in time. These 
are the important values to the thermal designer. The values are independent 
of the spacecraft and dependent only on the character of the earth-atmosphere 
system. The time average effective values have a probability distribution that 
is dependent on the time increment. A statistical study of current satellite 
data can give this probability distribution. With this distribution, the thermal 
designer can define confidence levels on predicted temperature ranges which 
are compatible with engineering models for use in design, failure probabilities, 
and spacecraft cost estimates. 
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TIME (t) 


LIMITATIONS OF p e (t) (HENCE INCLUDES p e (t) = p‘ e (At = o)), 
FOR 3 a LIMIT WORST CASE: 



< 

MAXIMUM 


C 0 + C,e 


-kAh 


( 2 ) 


P e (At) 



P e (t) dt/At 


< A., + a-) 

^ A 2 “ ■ ^2 


e K 1 At 
e K 2 At 


(3) p e ( At=o) =P e ( t) 


< .793 
>.080 


Figure 29. The limitations on the TAEA and an example of using 
these limitations to form a time history of the TAEA. 
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TABLE 9. TIME AVERAGE EFFECTIVE ALBEDO 


Sphere: At = 30 min Season = July Inclination = 35 deg 

Altitude (km) 

1 

2 

3 

500 




750 




1000 


+ 0. 19 
0-30 „ , , 
-0. 115 



The parameters for the thermal design engineer in the area of thermal 
and albedo radiation are summarized below: 

1. Probability distribution functions for the time average effective 
albedo and thermal radiation and the time rate of change distributions. 

2. The global albedo and thermal radiation averages. 

3. The functional dependence of the probability distribution functions of 
the variation of height, time increment, and solar zenith angle. 

4. The correlation function between the thermal infrared and albedo 
radiation. 

5. The separation of the radiation into specular and diffuse compo- 
nents, and surface and atmospheric components. 

Before the particular parameter studies can be undertaken, the probability 
distribution function for the time average effective values must be determined 
in a meaningful manner. It is proposed that these functions are highly desirable 
and necessary for future spacecraft thermal design studies. 


56 



REFERENCES 


1. Smith, R. E.: Space and Planetary Environment Criteria Guidelines for 
use in Space Vehicle Development. NASA TM X-64627, 1971 Revision, 
Nov. 15, 1971, pp. 1-11, 2-30 through 2-33. 

2. Bartman, F. L.: The Reflectance and Scattering of Solar Radiation by 
the Earth. University of Michigan, NASA CR-83954, 1967. 

3. Meyer-Arendt, R.: Radiometry and Photometry: Units and Conversion 
Factors. Applied Optics, vol. 7, 1968, pp. 2081-2084. 

4. Christensen, E. M. ; Radiation Geometric Factors Between Planets and 
Space Vehicles. AF 18( 600) - 1775, General Dynamics, 1962. 

5. Wolfe, W. L. , Ed.: Handbook of Infrared Technology. U. S. Government 
Printing Office, 1965. 

6. Earth Albedo and Emitted Radiation, NASA Space Vehicle Design 
Criteria (Environment). NASA SP-8067, July 1971. 

7. Snoddy, W. C. : Irradiation Above the Atmosphere Due to Rayleigh 
Scattering and Diffuse Terrestrial Reflections. AIAA Paper No. 65-666, 
AIAA Thermophysics Specialist Conference, Monterey, California, 

Sept. 13-15, 1965. 

8. Jones, R. W. : Factors Affecting the Near-Earth Thermal Environment. 
SMSD-SSL-1129, Teledyne Brown Engineering, Huntsville, Ala., 

Apr. 1970. 

9. Barnes, J. C. and Chang, D.: Accurate Cloud Cover Determination 
and Its Effects on Albedo Computations. Report No. 9G53-12, Allied 
Research Associates, Inc., Geophysics and Aerospace Division, 

Oct. 1968. 

10. Whitehill, L. P. : Survey on Earth Albedo. Tech. Rept. No. 1966-53, 
Lincoln Laboratory, MIT, Oct. 11, 1966. 

11. Pearson, B.D., Jr. and Neal, C. B.: Albedo and Earth-Radiation 
Measurements from OSO-II. Paper 67-330, AIAA Thermophysics 
Specialist Conference, New Orleans, La., Apr. 17-20, 1967. 


57 



REFERENCES (Concluded) 


12. Linton, R.C.: Earth Albedo Studies Using Pegasus Thermal Data. Prog, 
in Astron. and Aeron. , vol. 20, 1967, pp. 475-490. 

13. Raschke, E.: The Radiation Balance of the Earth-Atmosphere System 
from Radiation Measurements of the Nimbus II Meteorological Satellite. 
NASA TN D-4589, July 1968. 

14. Kornfield, J. and ITasler, A.F. : A Photographic Summary of the 
Earth' s Cloud Cover for the Year 1967. J. of App. Meteor. , vol. 8, 
no. 4, Aug. 1969, pp. 387-700. 

15. Stevenson, J. A.: Radiation Heat Transfer Analysis for Space Vehicles. 
ASD Technical Report 61-119, Space and Information Systems Division, 
North American Aviation, Inc. , Wright-Patterson Air Force Base, 

Ohio, Dec. 1961. 

16. Goldman, D. T. and Singer, S. F.: Studies of a Minimum Orbital 
Unmanned Satellite of the Earth (MOUSE), Part III Radiation and Tem- 
perature. Astronaut. Acta., vol. Ill, no. 2, 1957, pp. 110-129. 


58 



APPROVAL 


THE DEFINITION AND SPECIFICATION OF THE 
NEAR-EARTH ENVIRONMENTAL CRITERIA 
FOR SPACECRAFT THERMAL DESIGN 


By G. A. Gary and P. D. Craven 

The Information in this report has been reviewed for security classi- 
fication. Review of any information concerning Department of Defense or 
Atomic Energy Commission programs has been made by die MSFC Security 
Classification Officer. This report, in its entirety, has been determined to 
be unclassified. 

This document has also been reviewed and approved for technical 
accuracy. 



William c. snoddy 

Chief, Astronomy and So Ufa State Physics Division 



Director, Space Sciences Laboratory 


☆ U.S. GOVERMENT PRINTING OFFICE 1974 — 640-445 / 78 REGION NO. 4 


59 



DISTRIBUTION 


INTERNAL 
DA 01 

Dr. Lucas 
DS30 

Dr. Stuhlinger 
ES01 

Dr. Lundquist 
Mrs. Tingle (5) 

ES3 1 

Dr. Naumann 
ES2 1 

Dr. Decher 
ES11 

Mr. Snoddy 
Mr. Jones 

ES 12 

Mr. Arnett 
ES 14 

Mr. Miller 
Dr. Gary ( 15) 
Mr. Craven ( 15) 

ES4 1 

Mr. Vaughan 


CC01 

Mr. Wofford 
AS61 (2) 

AS61-L (8) 

AT01 (6) 

EXTERNAL 

Scientific and Technical Information 
Facility (25) 

P. O. Box 33 

College Park, Maryland 20740 

Attn: NASA Representative (S-AK/RKT) 

Joel M. Dash 
P. O. Box 8555 
General Electric Company 
Valley Forge Space Center 
Philadelphia, PA 10101 

Thomas H. Vonder Haar 
Dept, of Atmospheric Science 
Colorado State University 
Fort Collins, CO 80521 


R. L. Mancuso 
Stanford Research Institute 
333 Ravenswood Ave. 

Menlo Park, CA 94025 

Jaun K. Lovin 

Lockheed Missiles % Space Co. , Inc. 
Huntsville Research & Engineering Center 
4800 Bradford Drive 
Huntsville, AL 35807 

James C. Barnes 

Environmental Research & Technology, Inc. 
429 Marrett Road 
Lexington, MA 02173 

F. L. Bartman 

College of Engineering 

Dept, of Aerospace Engineering 

High Altitude Engineering Laboratory 

University of Michigan 

Ann Arbor, MI 48104 



